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ABSTRACT 

Based on the hard-scattering factorization which decomposes the diffractive 
structure function into a pomeron flux and a pomeron structure function, we study 
the single diffractive A+ production in polarized pp scattering: p + p p + A+ + X, 
which will be observed at forthcoming RHIC experiment. By analyzing the cross 
section and correlation of the spin polarization between the initial proton and 
produced A^, we found that the process might be effective for testing both hard- 
scattering factorization and models of the polarized gluon distribution in the pro- 
ton. 
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As is well known, the diffractive interaction which is characterized by a large 
rapidity gap event is described by pomeron exchange in the Regge theory 
Although the pomeron carrying the vacuum quantum number plays a crucial role 
in diffractive scattering, nature of the pomeron is still mysterious in the framework 
of quantum chromodynamics (QCD). Ingelman and Schlein suggested that the 
single pomeron exchange in proton-proton interaction can be probed as a hard 
scattering process between a hard parton in the pomeron emitted from a proton 
and a parton in another proton The evidence of the hard partonic structure 
of the pomeron was first reported by the UA8 Collaboration at CERN pp collider 
at y/s = 630 GeV by observing diffractive dijet production 0]. This result has 
been confirmed by further experiments ^, ^. In addition, recent experiment 
found the hard-gluon fraction in the pomeron to be 0.7 ± 0.2 [Q, which means 
the pomeron being an almost gluonic object. Those results suggest us to apply 
the factorization theorem to diffractive hard process as well as to usual inclusive 
processes. However, Collins proved that the factorization theorem for diffractive 
deep inelastic scattering (DDIS) is not expected to be applicable to hadron-hadron 
collisions, though this theorem works for the lepton and direct photon induced 
hard DDIS 0]. In fact, the predicted cross sections for hadron-hadron collisions are 
several times larger than the experimental data, and hence the hard factorization 
for DDIS does not work well for the hadron-hadron collision [§, ^j. In order to 
overcome such difficulty, Goulianos proposed a phenomenological model in which 
the structure of the pomeron is derived from the structure of the parent hadron ^j; 
he introduced a renormalized pomeron flux factor which was given by renormalizing 
the standard pomeron flux carried by a proton to be unity (renormalized pomeron 
flux model). Although the prediction by this model seems to be in agreement 
with experimental data [§, ^ , further tests of this model are necessary for various 
processes |p!0[ . 

On the other hand, the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory (BNL) will start soon. One of the important purposes of 
RHIC experiment is to extract information about the polarized gluon distribution 
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in the proton. That information is a key to understand the proton spin puzzle 
which is still one of the most challenging current topics in the nuclear and particle 
physics [|rT|. So far there are many models of the polarized parton distribution in 
the proton, which are extracted from a fit to the data on gi{x, Q^) [0, |T^. Those 
models can excellently reproduce experimental data on the polarized structure 
function of nucleons. However the behavior of the polarized gluon distributions is 
quite different among these models. In other words, the data on polarized structure 
functions of nucleon and deuteron alone are not enough to distinguish the model 
of gluon distribution functions. Although various processes have been proposed so 
far to test the models of the polarized gluon, knowledge of it is still poor. 

In this work, we propose another diffractive semi-inclusive process: p + p —* 
p + +X ( left side of Fig. [l| ), which can be a test of the renormalized pomeron 
flux model and also the polarized gluon distribution function. In this process, 
is dominantly produced via gluon-gluon fusion at the lowest order as shown in the 
right side of Fig. , where one of the 2 gluons is originated from the pomeron. 
Moreover, since is composed of a heavy quark c and antisymmetrically combined 
light u and d quarks, the spin of A+ is expected to be carried by the c quark 
and thus, the production of A+ in polarized proton-proton collisions gives us an 
interesting information about polarized gluons in the initial proton |n, |T^, rB|. To 



test the hard-scattering factorization and polarized gluon distribution model for 
this attractive process; 

P (Pj +P (Pj ^ P + A+ {p ) + X, (1) 
whose lowest order subprocess is0 

g(Pa) + g(pj c{pc) + c{pc), (2) 

where pi {i = A,B,a,b,c,c,A'^) in the parentheses denote the four-momentum 
of each particle and the over-arrow means that an initial gluon with momentum 



Since the charm quark content is extremely tiny and furthermore the pomeron is dominantly 
composed of gluons the gluon-gluon fusion is the dominant process for charm quark pair 
production. 

A subprocess in this model is a hard gluon-gluon scattering in the pomeron-proton system, 
as shown in Fig. |[ 
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Figure 1: The diagram for p + p ^ p + + X at the lowest oder in the framework 
of hard-scattering factorization (left side) and its subprocess diagrams (right side). 
In the left side, the particles with arrows are polarized. IP denotes the pomeron. 

Pb and a produced c quark with momentum pc are polarized, we introduce two 
useful observables; one is the dAa/dp^ defined in Eq.(^, which we call hereafter 
the polarized differential cross section, and the other is the All defined in Eq.(D, 
which we call the spin correlation asymmetry; 

dAa _ da{++) - da{-\ — ) + da{ ) - da{ — h) 



dpj, dprj, 

[da{++) - da{^ — ) + da{ ) - da{ — h)] /dp^ 

-) + dcx{-+)] /dp^ 



(3) 
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da/dp^ ' 

where da{^ — )/dp^, for example, denotes the spin-dependent differential cross sec- 
tion with positive helicity of the target proton and negative helicity of the produced 

Let us consider the process in the proton-proton center-of-mass frame. In this 
frame, we can take four-momenta of Pi as follows; 



Pa,b 

Pa+ = iE^+,p^,p 
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where the first, second and third components in parentheses are the energy, the 
longitudinal momentum and the transverse momentum, respectively, rrii is the 
mass of the i-particle. The polarized differential cross section, dAa/dp^, can be 
calculated as follows; 

dt dC, dQ / dx^ / dx^ 



dp^ J-l Jgmin Je™'" ^X^"" Jx'^'''' 

X f^%{^, t)f,/jp{xa, Q')AGg/,ix,, Q') 
x^^ADx+/,(^), (6) 

where t is the square of the four-momentum transfer of the proton which emits the 
pomeron. J' is the Jacobian which transforms the variables z and i into Q and p^. 
/g/jp(x^, Q"^) and /^^(^, t) are the hard gluon distribution function in the pomeron 
and the renormalized pomeron flux in the proton, respectively. The renormalized 
pomeron flux which was proposed by Goulianos to predict the observed single 
diffractive cross section f^, is defined by 

f^%iU)^Dfjp/,{U), (7) 

where fip/p{C,,t) and D are the standard pomeron flux and renormalization factor, 
respectively. fr/p{^, t) is given by 

/^/,(e,t) = ire'-2"WF2(t), (8) 

with the parameters which are chosen as [§] 

K = 0.73 GeV^ a{t) = 1 + 0.115 + 0.26[GeV-^]t, 

The renormalization factor D is defined as P] 

D = min(l,l) (9) 



with 



N = dd f]P/p{U)dt (10) 

JA/2/s Jt=Q 
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where Mq is 1.5 GeV^ being the effective diffractive threshold and the upper hmit 
for the ^ integration is a coherence hmit. In addition, AGg/^{xb, Q"^) and AD^+y^(z) 
represent the polarized gluon distribution function in the proton and the polarized 
fragmentation function of the outgoing charm quark decaying into a polarized Aj!", 
respectively. 

Since the subprocesses considered here are the same with the ones discussed 
in Ref. [14], here we repeat some important formulas for reader's convenience. By 
using the kinematical variables in Eq.(^, the polarized differential cross section, 
dAa/dt, for the subprocess is calculated to be 
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(11) 

t\Ux J 8 [ s ' 

where s, ti and U\ are defined as s = (p^ + pj^, ti = (p^ — p^)^ — and 

U\ = {p^ — pj^ — respectively. The Jacobian J' is given by 

2s/3p^ cosec^G , , 

J = , (12) 

z{s — 2mp)Jmj^+ +p^cosec^9 

where s and z are defined as s = {p^ + PbY = ^ + ' respectively, with 

•^1 — s-2mj, ' — s-2mj, ' 

In order to estimate the asymmetry. All, we need the unpolarized cross sec- 
tion which can be obtained by replacing polarized functions, AGg/^, dAa/di and 
AD^+^^by unpolarized functions, Gg/p, da/di and Da+/c' respectively, in Eq. (|]). 
The explicit formula of the unpolarized differential cross section for this subpro- 



cess was given by Babcock et al. [17|. As for the gluon distributions, we take the 
A AC [|13[ and the GRSVOl [1^ parameterization models for the polarized distribu- 
tion and the GRV98 |1^ model for the unpolarized one. Here, we set as (2mc)^ 
for each models. As for the unpolarized fragmentation function, we use Peterson 
fragmentation function |2^, D^^^+{z). However, unfortunately at present there 



are no established polarized fragmentation function AD^_^j^+ because of lack of 



experimental data. Thus, by analogy with the study on A polarization 
take the following ansatz: 



we 



- fi - 



where C'^^a+ scale- independent spin transfer coefficient. Here we apply the 
analysis on A production to production and choose the following two typical 
models: 

(i) C^__^yv+ ~ 1 (non-relativistic quark model), 

(ii) C^^A+ — ^ (J^t fragmentation model |]22[| ). 



For the hard gluon distribution function of the pomeron, we use R p3 



a;/g/7p(x, Q') = /g6x(l - x), /g = 0.7 ± 0.2, (13) 

where /g is the hard gluon fraction in the pomeron. To determine the value of the 
renormalization factor D given by Eq.(7), we estimated the value of defined by 
Eq.(|10D to be A^ = 3.4 (5.0) for ^ =200 (500) GeV using Eq.(|) 0. 

For numerical calculation, we use, as input parameters, rric = 1.25 GeV, nip = 
0.938 GeV and m^^ = 2.285 GeV 0. In numerically integrating Eq. (|), the 
minimum values of £, x and are given by f ™™ = x™™ = ttt^ — r and x™™ = 

^x^-xi ' respectively. In addition, to reduce possible non-pomeron contribution, we 



set ^™^''=0.05 as usual ^, |T0|. Furthermore, in order to get rid of the produced 
A+ entering the beam pipe, we limit the integration region of O for produced A+ 
as I < O < Pj, region is kinematically constrained from the condition of the 
diffiactive production. 

We show the p^. distribution of \d{A)(T/dpj.\ and All in Fig. | for = 200 
GeV and in Fig. ^ for ^/s = 500 GeV, respectively. Notice that the absolute value 
of dAcr/dpj. is presented in these figures, because the negative value of d/S.a/dp^ 
cannot be depicted in the figure which has an ordinate with logarithmic scale. 
Actually, the value of d/S.a/dp^ is negative for p^, region larger than the value 
corresponding to the apparent sharp dip shown in Fig. ^. Information on the 
dip is expected to be useful for distinguishing the models of the polarized gluon 
distribution functions because each parameterization model has a dip at different 
Pj, as shown in the left panel of Figs. ||. At -y/i = 500 GeV, the dip is not seen 
as shown in the left panel of Fig. ^ because the dip is in the region smaller than 
Pj, = 1 GeV for the case of y/s = 500 GeV. Thus, the polarized differential cross 
section are actually negative in the kinematical region presented in Fig. ^. From 
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Figure 2: The unpolarized and polarized differential cross section (left panel) 
and the spin correlation asymmetry (right panel) as a function of at ■\/s = 200 
GeV. The solid line in the left panel represents the unpolarized differential cross 
section with the GRV98 model for the unpolarized gluon distribution. The long- 
dashed and dashed lines show the polarized differential cross section calculated 
with C^^j^+ = 1 and C'c^a+ ~ respectively, for AAC and GRSVOl models of 
polarized gluon distributions. The value of dAa/dp^, is negative for region larger 
than the value corresponding to the apparent sharp dip. The same combination of 
the models for polarized gluon distributions and polarized fragmentation functions 
is adopted in the right panel. 
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Figure 3: The same as in Fig. 0, but for y/s = 500 GeV. 
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the right panel of Figs. and it seems that the All is an effective observable 
to distinguish various models of polarized gluon distribution functions, even if 
there are uncertainties of the spin-dependent fragmentation function. Notice that 
if the renormalized pomeron flux model is not taken into account for the present 
process, the polarized differential cross section and unpolarized differential cross 
section become 3.4 (5.0) times larger than our calculation for ^/s = 200(500) 
GeV. Therefore, measurement of those cross sections can be a good test of the 
renormalized pomeron flux model, though the All is not useful for testing the 
renormalized flux model because the renormalized flux factors are canceled out 
between the numerator and the denominator of Eq.(^. 

Finally, some comments are in order, regarding uncertainty of the calculated 
result on the value of rric and the choice of Q^. We have examined the variation 
of All on rric for the region of 1.15GeV < rric < 1.35GeV and found that the 
result remain to be unchanged. We also found that the dependence is rather 
insensitive; we have examined the 3 cases, (i) = (2mc)^, (ii) = p"^ and (iii) 

= without much difference. 

In summary, we have calculated the polarized differential cross section, d^a /dp^,, 
and the spin correlation asymmetry. All-, for the single diffractive production 
in polarized pp reactions at ^/s = 200 GeV and y/s = 500 GeV, based on the 
hard-scattering factorization with the renormalized pomeron flux. We found that 
dAa/dp^ and All largely depend on the model of the polarized gluon distribution 
function. Therefore, the process looks promising for testing the models of polar- 
ized gluon distribution function. Moreover, it is expected that the measurement 
of d{A)a/dp^ is quite effective for testing the renormalized pomeron flux model. 
In order to get more reliable prediction, the next-to-leading order calculation and 
error estimation might be necessary, which will be given in the forthcoming pa- 
per. In addition, since our prediction somewhat depends on the model of polarized 
fragmentation function, further experimental and theoretical study on polarized 
fragmentation functions is necessary for getting more reliable predictions. 

Although the present calculation is confined in the leading order, the results 
are interesting and we hope our prediction will be tested in the forthcoming RHIC 
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experiment. 
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